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Abstract 



The Dirac equation is solved for a pseudoscalar Coulomb potential in a 
two-dimensional world. An infinite sequence of bounded solutions are ob- 
tained. These results are in sharp contrast with those ones obtained in 
3+1 dimensions where no bound-state solutions are found. Next the general 
two-dimensional problem for pseudoscalar power-law potentials is addressed 
consenting us to conclude that a nonsingular potential leads to bounded solu- 
tions. The behaviour of the upper and lower components of the Dirac spinor 
for a confining linear potential nonconserving- as well as conserving-parity, 
even if the potential is unbounded from below, is discussed in some detail. 



The Coulomb potential of a point electric charge in a 1+1 dimension, 
considered as the time component of a Lorentz vector, is linear and so it 
provides a constant electric field always pointing to, or from, the point charge. 
This problem is related to the confinement of fermions in the Schwinger and 
in the massive Schwinger models 0-0 and in the Thirring- Schwinger model 
0. It is frustrating that, due to the tunneling effect (Klein's paradox), there 
are no bound states for this kind of potential regardless of the strength of 
the potential 0-0 • The linear potential, considered as a Lorentz scalar, 
is also related to the quarkonium model in one-plus-one dimensions 0-0. 
Recently it was incorrectly concluded that even in this case there is solely one 
bound state 0. Later, the proper solutions for this last problem were found 
0-^1]. However, it is well known from the quarkonium phenomenology in 
the real 3+1 dimensional world that the best fit for meson spectroscopy is 
found for a convenient mixture of vector and scalar potentials put by hand 
in the equations (see, e.g., |12j). The mixed vector-scalar linear potential in 
1+1 dimensions was recently considered |13j . There it was found that there 
are analytical bound-state solutions on condition that the scalar component 
of the potential is of sufficient strength compared to the vector component 
(\Vs\ > \Vt\)- As a by-product, that approach also showed that there exist 
relativistic confining potentials providing no bound-state solutions in the 
nonrelativistic limit. Although the discussion was confined to the vector- 
scalar mixing, the inclusion of a pseudoscalar potential could also be allowed. 

In a recent paper, McKeon and Van Leeuwen ^3] considered a noncon- 
serving-parity pseudoscalar Coulomb (NCPPC) potential (V = A/r) in 3+1 
dimensions and concluded that there are no bounded solutions for the reason 
that the different parity eigenstates mix. Furthermore, they asserted that the 
absence of bound states in this system confuses the role of the n-meson in 
the binding of nucleons. Such an intriguing conclusion sets the stage for 
the analyses by other sorts of pseudoscalar potentials. A natural question 
to ask is if the absence of bounded solutions by a NCPPC potential is a 
characteristic feature of the four- dimensional world. With this in mind, we 
approach in the present paper the less realistic Dirac equation in one-plus-one 
dimensions with a NCPPC potential (V = mu)\x\). 

The two-dimensional Dirac equation can be obtained from the four-dimen- 
sional one with the mixture of spherically symmetric scalar, vector and 
anomalous magnetic-like (tensor) interactions. If we limit the fermion to 
move in the x-direction (p y = p z = 0) the four- dimensional Dirac equation 
decomposes into two equivalent two-dimensional equations with 2-component 
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spinors and 2x2 matrices [T^]. Then, there results that the scalar and vector 
interactions preserve their Lorentz structures whereas the anomalous mag- 
netic interaction turns out to be a pseudoscalar interaction. Furthermore, 
in the 1+1 world there is no angular momentum so that the spin is ab- 
sent. Therefore, the 1+1 dimensional Dirac equation allow us to explore 
the physical consequences of the negative-energy states in a mathematically 
simpler and more physically transparent way. The confinement of fermions 
by a pure conserving-parity pseudoscalar double-step potential jTHj and their 
scattering by a pure nonconserving-parity pseudoscalar step potential [T7] 
have already been analyzed in the literature providing the opportunity to 
find some quite interesting results. Indeed, the two-dimensional version of 
the anomalous magnetic-like interaction linear in the radial coordinate, chris- 
tened by Moshinsky and Szczepaniak [18] as Dirac oscillator, has also received 
attention. Nogami and Toyama Toyama et al. j2U] and Toyama and 
Nogami [21] studied the behaviour of wave packets under the influence of 
that conserving-parity potential whereas Szmytkowski and Gruchowski [22] 
proved the completeness of the eigenf unctions. More recently Pacheco et al. 
[23] studied some thermodynamics properties of the 1+1 dimensional Dirac 
oscillator. 

Let us begin by presenting the Dirac equation in 1+1 dimensions. In 
the presence of a time- independent potential the 1+1 dimensional time- 
independent Dirac equation for a fermion of rest mass m reads 

H^ = E^ (1) 

H = cap + (3mc 2 + V (2) 

where E is the energy of the fermion, c is the velocity of light and p is the 
momentum operator, a and (3 are Hermitian square matrices satisfying the 
relations a 2 = (3 2 = 1, {a, (3} = 0. From the last two relations it steams that 
both a and (3 are traceless and have eigenvalues equal to ±1, so that one 
can conclude that a and (3 are even-dimensional matrices. One can choose 
the 2x2 Pauli matrices satisfying the same algebra as a and (3, resulting in a 
2-component spinor The positive definite function \^\ 2 = ty'ty, satisfying 
a continuity equation, is interpreted as a probability position density and its 
norm is a constant of motion. This interpretation is completely satisfactory 
for single-particle states [21]. We use a = <j\ and (3 = <t 3 . For the potential 
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matrix we consider 



V = lV t + (3V s + aV e + (3fV p (3) 

where 1 stands for the 2x2 identity matrix and (3^ = This is the most 
general combination of Lorentz structures for the potential matrix because 
there are only four linearly independent 2x2 matrices. The subscripts for the 
terms of potential denote their properties under a Lorentz transformation: t 
and e for the time and space components of the 2-vector potential, s and p 
for the scalar and pseudoscalar terms, respectively. It is worth to note that 
the Dirac equation is covariant under x — > —x if V e (x) and V p (x) change 
sign whereas V t (x) and V s (x) remain the same. This is because the parity 
operator P = exp(ie)P C3, where e is a constant phase and Po changes x 
into —x, changes sign of a and Pj 5 but not of 1 and (5. 
Defining the spinor ip as 

V> = exp ( % -k j (4) 



where 



A(x) = ^ dx'^p- (5) 
the space component of the vector potential is gauged away 

p + — J * = exp ( -^A J pip (6) 



c J \fo 

so that the time-independent Dirac equation can be rewritten as follows: 

Hip = Eip (7) 

H = o x cp + a 2 V p + a 3 (mc 2 + V s ) + lV t (8) 

showing that the space component of a vector potential only contributes to 
change the spinors by a local phase factor. 

Provided that the spinor is written in terms of the upper and the lower 
components 

*-(*) m 

the Dirac equation decomposes into : 
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(y t - E + V s + mc 2 ) 4>(x) = ihcx (x) + iVpx(x) 

(10) 

(Vt — E — V s — mc 2 ) x( x ) — ihc(j)'(x) — iV p (j)(x) 

where the prime denotes differentiation with respect to x. In terms of <p 
and x the spinor is normalized as f_°° dx (\<j)\ 2 + \x\ 2 ) — 1 ; so that <fi and 
X are square integrable functions. It is clear from the pair of coupled first- 
order differential equations (|10p that both and x must be discontinuous 
wherever the potential undergoes an infinite jump and have opposite parities 
if the Dirac equation is covariant under x — > —x. In the nonrelativistic 
approximation (potential energies small compared to the rest mass) Eq. (fTUj) 
loses all the matrix structure and becomes 



/ h 2 d 2 | ^ | y | v 2 | nv; \ (j) _r E mc 2\ ^ (12) 

\ 2m dx 2 s 2mc 2 2mc J 

Eq. (|TT|) shows that x if °f order t>/c << 1 relative to and Eq. (fT^j) 
shows that obeys the Schrodinger equation without distinguishing the con- 
tributions of vector and scalar potentials (at this point the author digress to 
make his apologies for mentioning in former papers (^1-^21) that the pseu- 
doscalar potential does not present any contributions in the nonrelativistic 
limit). It is remarkable that the Dirac equation with a nonvector potential, 
or a vector potential contaminated with some scalar or pseudoscalar cou- 
pling, is not invariant under V — > V + const., this is so because only the 
vector potential couples to the positive-energies in the same way it couples 
to the negative-ones, whereas nonvector contaminants couple to the mass of 
the fermion. Therefore, if there is any nonvector coupling the absolute val- 
ues of the energy will have physical significance and the freedom to choose a 
zero-energy will be lost. This last statement remains truthfully in the non- 
relativistic limit if one considers that such a contaminant is a pseudoscalar 
potential. 

Now, let us choose Vt — V s — and the intrinsically relativistic NCPPC 
potential V p = mujc\x\, where u is a real parameter. Although uj < gives 
rise to a potential unbounded from below, the possibility of such a sort of po- 
tential to bind fermions is already noticeable in the nonrelativistic limit of the 
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theory (see Eq. (|T2|l ). where the pseudoscalar linear potential, whether it is 
a conserving- or a nonconserving-parity potential, manifests itself effectively 
quadratic potential. Defining 



2m|u;| 

-li* (13) 

ff 2 - m 2 c 4 1 + g(cu) 
V ~ 2h\uj\mc 2 2 ^ 

where s(a>) stands for the sign function, the Dirac equation turns into 
the Schrodinger-like differential equations 



[v + 1/2 + a{u)]<f>, r]<0 



drf 4 



(15) 



d 2 x , rj 2 

+ tX 



[v + 1/2 + s{u)\ X , V>0 



d7]2 4 \(u+l/2) X , V<0 

The second-order differential equations (fl3j) have the form 



z 2 



V"(z) - [j + a J y(z) = 0, (16) 

whose solution is a parabolic cylinder function |25j. The solutions D^ a ^y 2 (z) 
and -D- a -i/2( — z) are linearly independent unless n = —a — 1/2 is a nonneg- 
ative integer. In that special circumstance D n (z) has the peculiar property 
that D n (—z) = (—l) n D n (z), and it is proportional to exp (— z 2 /A) H n (z/ y/2), 
where H n (z) is a Hermite polynomial. The solutions of (|T5|) do not exhibit 
this parity property so that we should not expect nonnegative integer values 
for v. The physically acceptable solutions for bound states must vanish in 
the asymptotic region \r}\ — ► oo and are expressed as 

= Q{- 7] )G^D v ^ ) {-^ + Q{+ n )G^D v {T ] ) 

(17) 

X = ei-^D^D^-^+ei+rfiD^D^M 
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where and are normalization constants and 9{rj) is the Heaviside 
function. Substituting the solutions (|T7j) into the Dirac equation (fTU|) and 
making use of the recurrence formulas 



—D v {z)-~D v {z) + D u+1 {z) 



^-D v (z) + Z -D u (z)-uD v ^(z) 



one has result 



(18) 



£>(+) 



£>(-) 



E + mc 2 
— mc 2 



x 



E + mc 2 
E — mc 2 



x 



- (z/ + 1) , > 

u; < 



■1/(^+1) 



> 

< 



The continuity of the wavefunctions (JT7J) at ?7 = furnishes 



(19) 



c(+) c(-) 



D(+) DH 



^+s(a;) 



(0) 



n 2 



A,(o) 



Together, (|T9*j) and lead to the quantization condition 



(20) 



A,+[i+sM]/2 (0) ± y/± + + s{to))/2]D uAl _ s{u))]/2 (0) = (21) 

This last result combined with ()14|) shows that the use of the minus sign 
inside the radical demands that — p < u < — lforu;>0, and — p < u < for 
uo < 0, here p stands for 1 + m 2 c 4 /2fi\uj\mc 2 . If, on the other side, one uses 
the plus sign then v > — 1 for uo > 0, and v > for uo < 0. 

The numerical computation of ()21)1 . which can be done easily with a 
symbolic algebra program, reveals no solutions for v < — 1. On the other 
side, for v > — 1 an infinite sequence of allowed values of v are found. The 
twenty lowest states for uo > are given in Table 1. The allowed values of 
v for uo < are those ones for oo > added by the unit. By inspection of 
Table 1 one sees that for v — > oo their values show a tendency to half-integer 
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numbers. The energy levels are obtained by inserting those allowed values of 
v in (ITU : 



E = ±^m 2 c i + 2h\ lu\ mc 2 {u + [1 + s(uS)\ /2} (22) 

One should realize that the energy levels are symmetrical about E = 0. It 
means that the potential couples to the positive-energy component of the 
spinor in the same way it couples to the negative-energy component. In 
other words, this sort of potential couples to the mass of the fermion instead 
of its charge so that there is no atmosphere for the production of particle- 
antiparticle pairs. No matter the intensity of the coupling parameter (uj), the 
positive- and the negative-energy solutions never meet. There is always an 
energy gap greater or equal to 2mc 2 , thus there is no room for transitions from 
positive- to negative-energy solutions. This all means that Klein's paradox 
does not come to the scenario. 

It is noticeable in ()17|) a somewhat left-right symmetry involving and 
X, such a symmetry is not exact inasmuch as C^' ^ D^ T \ The upper and 
lower components of the spinor have the same number of zeros, or nodes, and 
as a consequence of such a quasi-symmetry the zeros of and % exhibit, of 
course, an exact left-right symmetry. In whatever manner Eq. (|19 p -([20 p . and 
(|2*Tj) of course, are invariant under the change E — > —E with the proviso that 

the normalization constants transforms as |C^| < ► |-D^|, hence one can 

conclude that \<fi(±r)) \ < — > \x(Tv)\ i n suc h a wa Y that ^(±77) | — > lip^rj)]. 
If one recalls that the allowed values of v for uj > must be replaced by v + 1 
for uj < 0, a quasi-symmetry |0(±??)| — > |0(=p?)| and |x(±^)| — > |x(T^)| can 
also be observed in ()17|) under the transformation u — > —uj. This symmetry 
is not exact because neither |C^| — > |C^ T ^| nor — * Figures 1-4 

illustrate the behaviour of |0| 2 , \x\ 2 and \ip\ 2 = |0| 2 + |x| 2 for the ground-state 
contemplating all the possibilities of signs of E and uj. Comparison of these 
Figures shows that |?/>| tends to concentrate at the left (right) region when E 
and uj have equal (different) signs. It is also evident that |0| is larger (smaller) 
than |x| when E > (E < 0). Note that \<f>\ always tends to concentrate 
at the left when uj > (uj < 0). Figures 5 and Figure 6 illustrate the same 
things as Figure 1 for the first- and second-excited-states. Note from all these 
Figures that the probability position density has no nodes and the number 
of nodes of 0, which is equal to the one of \i increases with v. Furthermore, 
|0| is larger (smaller) than \x\ when E > (E < 0), independently of the 
sign of uj, as might be expected. 
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In addition to their intrinsic importance, the above conclusions renders 
a contrast to the result found in [14. There are bound-state solutions for 
fermions interacting by a pseudoscalar Coulomb potential in 1+1 dimensions, 
to tell the truth there is confinement, notwithstanding the spinor is not an 
eigenfunction of the parity operator. Therefore, the quadri-dimensional ver- 
sion of this problem requires clarification for the absence of bounded solu- 
tions. One might ponder that the underlying reason is the way the spinors 
are affected by the behavior of the potentials at the origin as well as at in- 
finity because this is the radical difference between the potentials in those 
two dissimilar worlds. In order to clarify this point let us consider (fTUj) for 
a general pseudoscalar potential. Those equations involve the coupling be- 
tween the upper and lower components of the spinor. This coupling can be 
formally eliminated when these equations are written as second-order differ- 
ential equations: 

H 2 d 2 $ 

' ' +\$ f * = E eff * (23) 



2m dx 2 eff 
where <3? refers to or x and 

E 2 — m 2 c 4 



E eff 



2mc 2 



(24) 



y± JL ± ^L +for ^ 

eff 2mc 2 2mc - for \ 



these last results show that the solution for this class of problem consists 
in searching for bounded solutions for two Schrodinger equations. It should 
not be forgotten, though, that the equations for (p or \ are n °t indeed in- 
dependent because the effective eigenvalue, E e ff, appears in both equations. 
Therefore, one has to search for bound-state solutions for V^ji with a com- 
mon eigenvalue. Now let us consider a nonconserving-parity pseudoscalar 
potential in the form V = fi\x\ s , then the effective potential becomes 

V± = J^-\ X \ 2S ± s ( x )M| x |*-i (25) 

Firstly let us consider the case \x > 0. When 5 > the effective potential 
goes to infinity as \x\ — > oo and it is finite at the origin for 5 > 1 whereas for 
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< 5 < 1 it is has a singularity given by ±s(x) fiS /2mc\x\ 1 s , implying for 
x > in a potential-well structure for Vjjj- and an attractive potential less 
singular than — l/8mx 2 for V e jj . The roles of V^j- and are changed for 
x < 0. Therefore, for 5 > the power potential leads to effective potentials 
fulfilling the conditions to furnish discrete spectra. On the other hand, when 
5 < the effective potential vanishes as \x\ — > oo and V^j (K//) ^ s always 
repulsive at the origin for x > (x < 0) in such a way that it is repulsive 
everywhere. Therefore, for 5 < the power-law potential does not lead to 
bound-state solutions. Note that these conclusions, either for 5 > or 5 < 0, 
are independent of the sign of /i on the condition that one has to change 
to V e jj if one changes the sign of ji. If one considers the conserving-parity 
potential V = jj,x s one obtains 

V ± - -^—x 25 ± ^x 5 - 1 (26) 
^ " 2mc^ ± 2mc X W 

and a similar analyzes also permits us to conclude that only when 5 > one 
could expect to find bounded solutions. 

Now then, we return to focus our attention on the linear potential. Both 
cases, if nonconserving- or conserving-parity, present potential-well struc- 
tures doing what is required to supply purely discrete spectra with infinite 
sequences of eigenvalues, in other words they are confining potentials, as it 
does every pseudoscalar power-law potential with 5 > 0. The nonconserving- 
parity potential V = mujc\x\ gives rise to the effective potential = 
mu 2 x 2 /2 ± s(x)fiw/2, an asymmetric potential with a discontinuity at x = 
equal to fkv. Now one can conceive why \ip\ and \<f>\ are more concentrated 
at the left region in the case E > and uo > 0. This behaviour is precisely 
what one would expect for such an asymmetric potential. On the other 
hand, the change u — > — to results in — > V^j and one should expect 
that l^l and \<j)\ be more concentrated at the right side. In addition, when 
the energy is sufficiently large (u large) compared to hw the discontinuity of 
the effective potential at the origin becomes unimportant and the fermion 
effectively undergo the influence of the average effective potential (mu 2 x 2 /2) 
maintaining their levels (in the sense of E e ff) nearly equally spaced. This 
interpretation is reinforced by observing Figures 1, 5 and 6, there one can 
see the propensity of to be closer to the origin for larger energies. As a 
matter of fact, a numerical calculation of the expectation value of i] (with 
m — u — c — h — 1) furnishes —0.505, —0.121 and —0.084 for the ground-, 
first-excited and second-excited-states, respectively. 
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The conserving-parity potential V = mujcx (the two-dimensional general- 
ized Dirac oscillator) is escorted by the effective potential V^jj = muj 2 x 2 /2 ± 
tku/2, which in addition to be continuous at the origin is even under x — > —x. 
In this case Eq. (JT5|) for 77 < must be disregarded and the equations valid 
for i] > must also be true for r) < 0, so the upper and lower components of 
the Dirac spinor are given by 

= e(- v )C^D„(- v ) +6(+ v )C^D u ( v ) 

(27) 

and the eigenvalues are even expressed by (J22J) . nevertheless the boundary 
conditions implies that C^~> = C^' and = D^ + \ with 



£>(+) 



2 



E + mc 2 f -(u + 1), uj>0 



E - mc 2 \ < 



(28) 



In addition, the boundary conditions implies, for uj > 0, that D u (0) = 
with allowed values v = 1, 3, 5, . . . and D' v (0) = with allowed values v = 
0,2,4, . . . In this case D v (—rf) = (—l) u D v (rj), and D v {rj) is proportional to 
exp (— r/ 2 /4) H u (r]/^/2), where H u (rj) is a Hermite polynomial (for uj < 
one has to change the allowed values of v by v + 1). The upper and lower 
components of the spinor have definite and opposite parities under rj — > —i] 
and, since the Hermite polynomial H v (rj) has v different zeros, the number 
of nodes of x is the number of nodes of plus (minus) one for uj > {uj < 0). 
This fact can be better understood by keeping in mind that V^jj = V^f±Huj. 
The energy levels (in the sense of E e ff) are always equally spaced because 4> 
as much as x are affected by continuous and symmetric harmonic potentials. 
The symmetry already remarked for the potential V = mujc\x\ translates 
into IC^I < — > l-D^I and \(fi(r])\ < — ► |x( r ?)|; under the change E — ► —E, in 
such a manner that is invariant. Also a quasi-symmetry is present under 
uj — > — uj: |0(77)| « — ► |x( ? ?)l- m this case the symmetry is not perfect because 
neither \C^\ -> nor -> All in all, \<j>\ is larger (smaller) 

than \x\ when E 1 > (E < 0), independently of the sign of uj. 

For short, in one-plus-one dimensions, we have succeed in searching bound- 
state solutions for the Coulomb potential as well as for the generalized Dirac 
oscillator. These problems are somehow similar as long as they are affected 
by effective harmonic potentials and they, in fact, possess the same formal 
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expression for the eigenvalues. Nonetheless, the quantum numbers and their 
related eigenfunctions differ because those different problems are subject to 
distinct boundary conditions 

As stated in the third paragraph of this work, the anomalous magnetic- 
like interaction in the four-dimensional world turns into a pseudoscalar in- 
teraction in the two-dimensional world. The anomalous magnetic interaction 
has the form —ifij3a.'V<j)(r), where \i is the anomalous magnetic moment in 
units of the Bohr magneton and <fi is the electric potential, i.e., the time 
component of a vector potential [21]. m one-plus-one dimensions the anoma- 
lous magnetic interaction turns into o"2/i0', then one might suppose that the 
pseudoscalar Coulomb potential is due to an electric potential proportional 
to s(x)x 2 /2. The oddness of the electric potential (under x — > —x) does 
not present problem to the confinement of a fermion because its effective 
mass, an x-dependent mass which always increases as the fermion goes away 
from the origin, depends on (</>') , an result independent of the sign of s(x) 

uni-mi 

A word should be said about the role of the 7r-meson field. The La- 
grangian density describing the pion-nucleon interaction C = —iXip^^jr is 
parity-invariant because the 7r-meson is a pseudoscalar field, i.e., 7r(r, t) — > 
— 7r(— r, t) under the parity transformation. Nonetheless, the interaction ma- 
trix term present in the Dirac equation as written by McKeon and Van 
Leeuwen [T3j, i(5^ b V(r), supposed to be due to the 7r-meson field is not 
parity-invariant due to the reason that the potential function, V(r), is parity- 
invariant but the potential matrix is not. This argument exposes the inad- 
equacy of the interaction term in the Dirac equation as proposed in |14j . 
Therefore, any conundrum about the role of the 7r-meson should be con- 
sistently presented by taking into account a quintessential parity-invariant 
potential. Furthermore, in order to correspond to the physical reality one 
should be aware that the massive 7r-meson field gives rise to a Yukawa po- 
tential instead of a Coulomb potential. 

The methodology of effective potentials prompted in this paper might 
be extended for the four-dimensional world. Such a methodology seems 
more difficult there because the upper and lower components of the Dirac 
spinor involve not only radial functions but also angular functions (the two- 
component spinor spherical harmonics). Probably, we have to deal with four 
coupled second-order differential equations. This task of carry the Dirac 
equation through equivalent Sturm-Liouville problems is under way and will 
be reported in due time. 
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Figure captions 

Figure 1 - \<p\ 2 (full thin line), \x\ 2 (dashed line) and \<p\ 2 + \x\ 2 (full 
thick line), corresponding to positive-gr ovjid-st&te energy for the potential 
V = mujc\x\ with uj > and m — \ou\ — c — h — 1. 

Figure 2 - The same as in Figure 1 with uj < 0. 

Figure 3 - The same as in Figure 1 corresponding to negatwe-ground-state 
energy {uj > 0). 

Figure 4 - The same as in Figure 1 corresponding to negatwe-ground-state 
energy with uj < 0. 

Figure 5 - The same as in Figure 1 corresponding to positive-first-excited- 
state energy (uj > 0). 

Figure 6 - The same as in Figure 1 corresponding to positive-second- 
excited-state energy (uj > 0). 
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Table 1: The lowest allowed values of is, for u > 0, such that Eq. (21) is 
satisfied. 



v for plus sign in front of the radical 


v for minus sign in front of the radical 


-0.654541 


0.5485571 


1.468573 


2.522295 


3.482395 


4.514350 


5.487785 


6.510563 


7.490650 


8.508353 


9.492428 


10.506906 


11.493638 


12.505886 


13.494514 


14.505129 


15.495179 


16.504543 


17.495700 


18.504078 
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